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-  ABSTRACT  - 

The  electronic  structure  and  geometric  characteristics  of 
small  Pd  clusters  (up  to  13  atoms)  are  studied  at  the  SCF/CI 
level,  with  major  focus  on  the  mutual  dependence  between  these 
properties.  To  this  end,  Jahn-Teller  effects  have  been 
considered,  and  the  magnetic  moments  analyzed  on  the  distorted, 
most  stable ,  structures . 

By  means  of  calculations  of  the  INDO  type,  the  energies  and 
normal  modes  involved  in  the  different  steps  of  the  distortions 
which  leads  to  diamagnetic  Pdi3  but  paramagnetic  Pd4  are 
described.  Comparisons  with  other  calculations  and  experimental 
data  are  included. 
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1-  INTRODUCTION  - 

The  potential  uses  of  metal  clusters  in  chemistry,  physics 
and  technology  have  lead  to  impressive  effort  towards  a  better 
understanding  of  their  properties.  Because  their  characterization 
represents  an  important  link  in  the  understanding  of  the 
fundamental  mechanisms  of  catalysis,  the  basic  properties 
(geometry,  bond  strength,  reactivity)  of  small  metal  aggregates 
have  become  the  subject  of  intense  theoretical  [1-13]  and 
experimental  study  [14-21]. 

It  is  well  known  that  the  electronic  properties  and 
structural  characteristics  of  metal  aggregates  are  different  from 
the  bulk  [20-22],  The  nature  of  these  differing  electronic 
properties  is  believed  to  play  an  important  role  in  the 
adsorption  and  reaction  of  small  molecules,  and  is  largely 
related  to  geometry  [6,16,23,24].  Metal  clusters  of  suitable  size 
are  able  to  reproduce  the  electronic  features  of  crystal  defects 
and  of  the  small  metal  particles  in  supported  structures  that  are 
responsible  for  catalytic  reactivity. 

Research  on  the  electronic  structure  and  geometry  of  small  metal 
clusters  may  lead  to  the  development  of  more  selective  and 
efficient  catalysts,  and  to  model  specific  reactions  through  the 
analysis  of  the  interactions  involved  in  the  adsorption  of  simple 
molecules . 

As  the  simplest  representation  of  the  reactive  site,  often 

metal  clusters  defined  by  2  to  100  atoms  are  studied,  with 

particular  emphasis  in  the  change  of  the  properties  as  a  function 

of  the  size  [3,4,6,10].  However,  as  has  been  previously 

described,  not  only  the  size,  but  also  the  geometric  and 

electronic  structure  are  important  to  reproduce  surface  defects 

or  metal  particles  in  supported  catalysts.  Besides,  the 

experimental  information  available  on  the  geometrical  structure  r 

of  metallic  clusters  of  the  order  of  10  atoms  is  indirect 

[15,20,21]  and,  in  most  cases, not  even  the  symmetries  are  known. 

Any  calculation  of  the  physical  properties  of  microclusters  has 

to  start,  therefore,  by  the  prediction  of  the  geometry  by  means 

of  the  minimization  of  the  total  energy.  On  the  other  hand,  the  / 

most  symmetrical  configurations  of  small  clusters  are  very  often  f  Cofles 
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unstable  due  to  the  Jahn-Teller  effect  [1,4,7,25],  The 
calculations  that  only  consider  highly  symmetric  geometries  often 
miss  the  true  equilibrium  geometry.  For  a  detailed  understanding 
of  the  structure-reactivity  relationships  in  clusters,  and  a 
successful  correlation  between  chemical  and  physical  properties, 
the  most  stable  structure  of  the  cluster  must  first  be 
determined.  However,  because  of  the  unquestionable  coupling 
between  electronic  and  geometric  structure,  clearly  shown  by  the 
importance  of  the  Jahn-Teller  distortions  [7],  both  are  to  be 
determined  simultaneously  when  an  interpretation  of  the 
reactivity  is  the  goal. 

This  mutual  dependence  follows  from  the  consideration  of  two 
coupled  effects: 

The  non-equivalent  occupation  of 
energetically  degenerate  orbitals  leads  to  a  splitting  of  the 
degeneracy,  by  means  of  a  distortion  of  the  geometry  to  a  lower 
symmetry  one.  The  structure  of  the  cluster  is  determined  by  Jahn- 
Teller  effects. 

Elemental  clusters,  on  the  other  hand,  have 
the  tendency  to  assume  highly  symmetric  structures,  as  a 
consequence  of  which  degeneracy  in  the  one-electron  levels  occur 
(effect  of  the  geometry  on  the  electronic  structure) .  Only  those 
systems  with  a  number  of  electrons  sufficient  to  completely  fill 
or  half-fill  the  sets  of  degenerate  one-  electron  levels  will  not 
Jahn-Teller  distort. 

This  electronic  structure-  geometry  interplay  opens  a 
challenging  field  in  the  study  of  the  characteristics  of  metal 
clusters,  when  properties  relevant  to  catalysis  (local  density 
charges,  magnetic  moments)  are  to  be  analyzed. 

With  the  aim  of  better  understanding  the  coupling  between 
the  electronic  configuration  and  the  geometry  of  the  transition 
metal  structures,  we  discuss,  in  this  paper,  the  electronic  and 
geometric  characteristics  of  Pd  clusters,  evaluated  at  the  Self 
Consistent  Field  -  Configuration  Interaction  (SCF/CI)  level,  in 
an  attempt  to  compare  them  with  those  previously  found  for  Rh 
[23-25]. 
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Different  cluster  sizes  can  be  analyzed  for  a  given 
transition  metal.  On  the  basis  of  previous  results  on  Rh  clusters 
[24],  we  have  chosen  a  13  atom  cluster  because  it  is  large  enough 
to  provide  a  good  description  of  the  behavior  of  the  tiny 
particles  that  are  present  in  supported  catalysts  and  small 
enough,  on  the  other  hand,  to  allow  accurate  calculations.  A  13 
atom  cluster  is  also  compatible  with  the  representation  of 
structures  of  different  possible  symmetries  [5,25],  and  13  is  one 
of  the  magic  numbers  of  atoms  frequently  observed  in  clusters 
of  various  kinds  [5,26,18,19].  In  addition,  going  from  the  bulk 
to  the  smaller  crystallites,  the  density  of  states  of  transition 
metal  13  atom  clusters  [27]  has  been  found  to  be  close  to  that 
of  the  bulk.  Finally  we  remark  that  often  properties  that  bear 
no  resemblance  to  the  crystalline  state  are  more  frequently 
found:  the  most  attractive  include  the  larger  stability  of  the 
five  fold  symmetry  structures  and  the  unusual  magnetic  properties 
of  these  systems. 

From  many  experimental  [18,19,27,28]  and  theoretical  [3- 
7,10,11]  investigations,  it  was  concluded  that  five  fold  symmetry 
is  the  natural  choice  of  small  microclusters  of  materials  that 
crystalize  in  fee  lattices.  Such  structures  are  based  on  a 
centered  13-atom  icosahedron,  representing  the  minimum  energy 
packing  configuration  in  mono-  and  bimetallic  clusters  [29-33]. 
The  planes  of  13  as  well  as  55-atom  icosahedral  structures 
resemble  close  packed  (111)  planes,  the  most  compact  and  most 
stable  single  crystal  structures  of  fee  lattices.  While  there  is 
evidence  that  a  55  atom  icosahedron  likely  represents  the  actual 
particle  size  in  supported  clusters,  a  13  atom  icosahedon  is  an 
obvious  simpler  target  that  keeps  the  same  local  properties  on 
the  reactive  centers  because  of  the  identical  local  environment 
on  each  of  the  sites.  This  provides  an  extra  justification  for 
the  choice  of  this  cluster  size  in  our  studies. 

According  to  the  previous  discussion  we  apply,  in  this  work, 
a  version  of  the  Intermediate  Neglect  of  Differential  Overlap 
(INDO)  model  [34,35]  at  the  Self  Consistent  Field-  Configuration 
Interaction  (SCF-CI)  level,  to  study  the  results  of  the  geometry- 
electronic  structure  coupling  on  the  magnetic  and  structural 
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characteristics  of  Pdi3  clusters.  Results  on  Pd4  and  on  the  dimer 
Pd2  are  also  reported  for  comparison  and  calibration. 


2-  COMPUTATIONAL  DETAILS  - 

Because  of  the  4d^®  (^Sg)  ground  state  of  the  Pd  atom,  open- 
and  closed-shell  electronic  configurations  are  compatible  with 
both  an  odd  or  an  even  number  atom  cluster.  The  stability  of  the 
different  multiplicities  (M)  for  a  given  cluster  size  has  been 
compared  after  a  complete  search  for  the  structure  of  minimum 
energy  on  the  potential  energy  surface,  starting  from  different 
initial  geometries  which  are  associated  with  the  different 
symmetries  compatible  with  the  number  of  atoms  that  define  the 
cluster. 

The  calculated  geometries  are  the  result  of  a  full 
optimization  of  the  coordinates  (interatomic  distances  and 
angles) ,  without  any  constraint  on  their  variation.  Optimization 
is  based  on  a  minimization  of  the  gradient  [36,37],  evaluated 
analytically,  using  the  BFGS  algorithm  to  update  the  inverse 
Hessian  matrix  in  successive  geometry  cycles.  Calculations  have 
been  made  at  the  Restricted  Hartree  Fock  (RHF)  and  Restricted 
Open-Shell  Hartree  Fock  (ROHF)  level.  Dealing  with  structures  of 
high  symmetry  and  a  large  number  of  electrons,  changes  in  the 
one-electron  distribution  during  the  SCF  cycles  may  lead  to  a 
non-equivalent  occupation  of  degenerate  orbitals  which,  by 
breaking  the  symmetry  in  the  electronic  distribution  may  result 
in  spurious  Jahn  Teller  distortions.  In  order  to  avoid  this 
effect,  we  start  the  calculations  using  Configuration-Average 
Hartree  Fock  (CAHF)  theory  [38],  with  an  average  M  for  the  number 
of  electrons  considered.  The  final  geometries  depend  somewhat  on 
the  number  of  open-shell  orbitals  and  electrons.  This  number  has 
been  chosen  to  be  compatible  with  the  size  of  the  open-shell  for 
the  different  average  M  (i.e.,  4  electrons  in  3  orbitals  for  M=3, 
but  in  6  orbitals  for  M=6)  .  For  each  of  the  different  geometries 
that  result  from  the  different  definition  of  the  open-shell,  the 
M  has  been  chosen  after  Cl  calculations,  using  the  orbitals  of 
the  CAHF  calculation  as  the  reference  for  a  Rumer  Cl  projection 
over  pure  spin  states  [39].  In  this  way,  M  and  geometry  are 
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simultaneously  optimized,  avoiding  the  possibility  of  getting 
trapped  in  spurious  local  minima  corresponding  to  a  stable 
geometry  for  a  given  electronic  configuration. 

For  the  13  atom  Pd  cluster,  the  stability  of  icosahedral 

(Ijj)  precrystalline  structures  has  been  compared  to  that  of  fee 

(Ojj)  and  hep  pieces  of  bulk.  In  order  to  keep  the  Ijj  symmetry,  it 

was  neccessary  to  start  with  an  average  of  130  electrons  in  78 

orbitals,  then  decreasing  the  orbital  space  to  66  using  the 

previous  set  of  eigenvectors  as  an  initial  guess.  These  vectors 

then  serve  as  starting  orbitals  for  a  calculation  with  4 

electrons  in  3  orbitals.  With  a  similar  strategy,  130  electrons 

* 

were  averaged  in  78  orbitals  to  start  the  SCF  calculations  for 
fee  (Oh)  and  hep  symmetries.  The  resulting  orbitals  were  used  for 
an  intermediate  calculation,  before  the  final  RHF,  with  130 
electrons  in  65  orbitals. 

In  the  case  of  Rh4,  geometry  optimizations  have  been  done 
for  closed-shell  and  an  average  of  4  electrons  in  3  orbitals,  4 
electrons  in  6  orbitals,  and  10  electrons  in  9  orbitals, 
corresponding,  in  this  way,  modelling  possible  multiplicities  (M) 
as  high  as  3,  5  or  7,  respectively.  With  this  strategy, 

distortions  in  the  geometry  during  the  SCF  cycles  are  avoided, 
and  geometries  optimized  for  different  electronic  configurations 
are  compared  after  a  Cl  that  uses  the  resulting  orbitals  that 
manta in  full  symmetry. 

It  is  clear  from  this  work  and  our  previous  work,  that 
correlation  energy  must  be  included  in  the  comparison  of  the 
energies  of  the  different  geometries  for  a  given  cluster  size. 
Final  energies,  as  well  as  electronic  parameters  and  molecular 
orbital  interactions,  are  the  result  of  multi-reference 
configuration  interaction  calculations  using  single  excitations, 
MRCIS. 

Two  parametrizations  of  the  INDO  theory  have  been  used  in 
these  studies;  one  for  geometry,  which  utilizes  two-center  two- 
electron  integrals  that  are  calculated  ab-initio;  and  one  to 
calculate  the  electronic  descriptors  and  compare  different  M  at 
fixed  geometries,  that  obtain  these  integrals  empirically  from 
atomic  spectroscopy  [35,  40-42].  This  latter  method  is 
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parametrized  directly  on  molecular  spectroscopy  at  the  Cl  singles 
(CIS)  level,  which  explains  our  use  of  MRCIS  for  these  clusters. 
The  resonance  integrals  P  are  chosen  according  to  formulae  that 
takes  into  account  different  electronegativities  [23],  and 
reproduces  the  available  experimental  geometries  and  spectroscopy 
for  the  dimer. 

More  details  on  the  method  are  given  elsewhere  [24,41]. 

3-  RESULTS  AND  DISCUSSION  - 

3.1-  Pdj  dimer  - 

The  ground  state  of  transition  metal  dimers  is  often 

* 

controversial,  because  both  the  symmetry  and  the  geometry  are 
largely  dependent  on  the  level  of  correlation  used  in  the 
calculations. 

Our  MRCIS-INDO  calculations  give  a  triplet  ground 

state  for  Pd2,  with  a  leading  configuration 

and  an  interatomic  distance  of  2.46  A.  The  lowest  exited  state 
is  a  singlet  (^Sg'*’)  ,  0.16  eV  higher  in  energy  than  the  first 
triplet,  with  a  Pd-Pd  bond  length  of  2.67  A.  The  next  triplet, 
^Ilg,  lies  0.44  eV  above  the  ground  state. 

The  Mulliken  population  analysis  (5s°-°®5p°-°®4d®-®^  in  the 
first  exited  singlet,  5s0*5i5p0*054jj9.44  triplet  ground 

state)  indicates  a  larger  contribution  of  the  5s  orbital  to  the 
stabilization  of  the  metal-metal  bond  in  the  triplet,  shortening 
the  bond  length  relative  to  the  singlet.  The  4d^°  Pd  atoms 
undergo  rehybridization  of  orbitals  in  order  to  form  the  5s 
bonding  combination  Og  that  defines  the  HOMO. 

Our  results  on  the  spectroscopy  and  the  geometry  of  the 
dimer  are  in  agreement  with  Density  Functional  [43,44]  and 
MCSCF/MRCISD  calculations  [45,46].  The  latter  has  demonstrated 
to  give  results  that  agree  well  with  experiment,  not  only  for  the 
dimer  Pd2  [45]  but  also  for  the  dimer  Pd3  [46]  .  There  is 
general  agreement  among  the  different  authors  in  the  well  studied 
case  of  Pd2,  in  contrast  to  the  situation  for  Rh2  [23].  While  a 
mixture  of  configurations  characterize  the  electronic  states  in 
Rh2  [23,47,48],  we  have  found  that  nearly  pure  configurations 
characterize  the  lowest  lying  states  of  Pd2.  There  is  also 
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agreement  in  the  larger  bond  length  in  Pd2  than  in  Rh2  (2.28  A 
[23]),  associated  with  a  more  efficient  5s-5s  bonding  in  the 
latter,  where  no  promotion  energy  is  needed  for  the  5s 
contribution  to  the  a  bond. 

3.2-  Pd^  structures  - 

Four  atom  clusters  are  most  frequently  studied  in 
tetrahedral  symmetry  as  models  for  a  (111)  structure  [44].  They 
are  usually  the  most  stable  geometries  for  fee  habits,  with  a 
large  coordination  number,  compatible  with  the  directionality  of 
the  d  orbitals.  We  have  compared  tetrahedral  and  square  planar 
structures,  and  found  the  former  1.77  eV  more  stable  than  the 
latter . 

The  final  geometry  of  the  tetrahedral  Pd4  is  strongly 
dependent  on  the  number  of  open-shells  included  in  the 
configuration  average,  as  shown  in  table  1,  and  this  is  to  be 
expected.  The  more  states  that  we  include  in  an  average  involving 
frontier  orbitals  and  electrons,  the  more  likely  we  are  to 
include  high  energy  states  in  the  average.  At  the  SCF  level,  RHF 
calculations  are  based  on  a  4d^°  configuration  of  the  atoms  and 
lead  to  longer  interatomic  distances  than  the  ROHF  procedure. 
These  calculations  imply,  through  the  definition  of  the  open- 
shells,  that  the  participation  of  the  5s  orbitals  is  important 
in  the  stabilization  of  the  bond.  Different  M  have  been  compared 
at  the  MRCIS  level  for  each  of  the  structures  in  table  1. 
Independent  of  the  specific  model  and  the  geometry,  the  triplet 
is  always  the  most  stable  state. 

In  table  1  the  geometry  of  the  triplet  is  determined  from 
the  indicated  calculation  using  the  INDO/1  model  with  parameters 
determined  for  geometry  as  described  in  the  Computational  Details 
section.  The  relative  energies  at  fixed  geometry  are  then 
determined  via  MRCIS  INDO/S  calculations.  For  this  reason 
energies  between  rows  cannot  be  compared. 

The  Mulliken  populations,  reported  for  the  triplets,  show 
the  larger  contribution  of  the  s  orbitals  to  the  shorter  bonds. 
The  s  participation  is,  in  turn,  dependent  on  the  open-shell 
definition  of  the  initial  configuration. 
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Lower  energies  are  always  found  for  an  average  of  4 
electrons  in  3  orbitals.  For  this  definition  of  the  open-  shell, 
several  interatomic  distances,  which  may  either  result  from 
geometry  optimizations  from  different  initial  guesses  or  single 
point  calculations,  have  been  compared.  As  one  of  the  single 
points,  the  interatomic  distance  characteristic  of  the  bulk  (2.75 
A  [49])  has  been  analyzed.  The  results  for  the  different  M  are 
reported  in  table  1.  These  calculations  indicate  a  very  flat 
potential  for  Pd4.  Within  the  accuracy  of  our  treatment,  however, 
no  change  in  the  interatomic  distance  with  the  change  in  M  can 
be  assessed. 

* 

The  most  stable  structure  (ground  state-GS)  corresponds  to 
a  triplet  ^T2  (Fig  1) ,  with  a  leading  configuration 
laj^2lt2®le^lti®2e^2ti®  2t2^2aj^3t2^3t^^ .  At  this  geometry  the 
singlet  ^A^  lies  0.40  eV  higher  in  energy  and  has  the  same 
electronic  assignment  as  the  GS. 

The  Mulliken  population  analysis  shows  a  smaller  population 
of  the  s  orbital  in  the  configuration  that  results  from  RHF 
calculations  and  this  leads,  in  turn,  to  longer  Pd-Pd  interatomic 
distances . 

Similar  calculations  for  the  square  planar  geometry  yield 
a  GS,  with  an  interatomic  distance  of  2.34  A. 

This  state  is  0.14  eV  lower  in  energy  than  the  first  exited 
state  we  calculate  for  this  structure,  which  is  of  symmetry . 

3.3-  Pd^3  structures  - 

After  a  full  optimization  of  the  geometry  in  the  way 
previously  described,  the  symmetry  is  calculated  to  be  most 
stable,  table  2.  The  interatomic  distance  of  2.63  A  among  the 
peripheral  surface  atoms  (Fig.  2)  represents  a  compromise  between 
the  value  in  the  diatomics  (2.46  A  [48]),  and  that  in  the  bulk 
(2.75  A  [49] ) . 

The  results  of  MRCIS  calculations  show  that  the  GS  is  a 
triplet,  ^T2g  (Fig.  3),  0.21  eV  lower  in  energy  than  the  first 
exited  singlet  ^Hg.  These  two  states  arise  mostly  from  the 
configuration  shown  in  Fig.  3.  The  next  triplet,  ^Gg,  is  0.28  eV 
above  the  GS,  followed  by  two  singlets,  ^Gg  and  ^Ag,  0.08  and 
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0.25  eV  higher  in  energy  respectively  .  The  singlet  ^Ag  also 
arises  from  the  configuration  shown  in  Fig.  3.  The  lowest  quintet 
^Ty,  lies  0.73  eV  above  the  GS  (table  2).  About  38  terms  lie 
within  1.0  eV  of  the  lowest  calculated  ^T2g. 

The  most  stable  structures  in  either  Oj^  or  03}^  symmetries 
belong  to  closed-shell  configurations  ^Ag  and  ^A^”,  respectively, 
and  are  about  1.7  eV  less  stable  than  the  ^T2g  GS  of  symmetry 
(table  2).  The  triplet  ^T2g  in  Ojj  symmetry  is  0.37  eV  higher  in 
energy  than  the  singlet  ^Ag,  and  the  quintet  of  ^T^y  symmetry  is 
1.15  eV  higher  in  energy.  In  the  hep  (D3jj)  structure,  the  triplet 
^Ai”  and  the  singlet  ^A^^”  are  almost  degenerate,  with  the  singlet 
0.06  eV  lower.  The  lowest  quintet  of  ^A^' symmetry  is  calculated 
to  lie  1.036  eV  above  the  GS. 

We  have  checked  the  ^Ag  of  Ojj  symmetry  and  the  ^A^"  of  D3jj 
symmetry  and  both  are  true  minimum  on  the  potential  energy 
surface. 

The  larger  stability  of  five  fold  symmetries  in  small 
transition  metal  clusters  has  been  both  experimentally  determined 
[15,20,21]  and  theoretically  derived  [5,7,11,24,50].  It  is  known 
from  experiment  that  five  fold  symmetry  characterizes  small  Co, 
Cu,  Ni  and  Mg  clusters  [15,20,21].  Calculations  are  more 
frequently  found  for  the  lighter  elements,  although  molecular 
dynamics  and  density  functional  theory,  applied  to  Ni23  and  CUj^3 
respectively  [3,5,11],  also  favor  the  Ij^  symmetry.  In  recent 
density  functional  calculations  (LSD) ,  Reddy,  Khanna  and  Dunlap 
[50]  find  Ijj  structures  more  stable  than  0^  for  Pd2^3,  Rh23  and 
RUj^3,  in  agreement  with  our  present  results  on  Pd  and  those 
previously  published  for  Rh  [25].  However,  they  have  not 
considered  any  decrease  in  the  symmetry,  and  the  large  magnetic 
moment  of  the  1^  clusters,  which  results  from  the  partial  filling 
of  the  d  bands  in  structures  of  high  symmetry,  is  the  main  topic 
of  their  discussion. 

The  small  transition  metal  clusters  are  models  of  catalytic 
active  sites,  which  are  usually  associated  with  surface  defects, 
characterized  by  a  local  electron  deficiency.  A  Mulliken 
population  analysis  (table  2)  shows  a  negative  charge  (-0.175  au) 
largely  concentrated  on  the  atom  in  the  center,  counterbalanced 
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by  a  small  positive  charge  (0.015  each  atom)  distributed  among 
the  peripheral  atoms.  Previous  calculations  on  the  interaction 
of  CO  with  Rh  clusters  [24]  have  demonstrated  larger  activity 
towards  CO  adsorption  on  electron  deficient  sites, 
characteristics  of  the  surface  atoms,  in  the  precrystalline 
structures.  The  study  of  the  electronic  characteristics  imply, 
thence,  practical  importance,  and  can  only  be  analyzed  after  a 
complete  study  of  the  structure  of  minimun  energy. 

3.4-  On  the  stability  of  the  highly  symmetric  structures  - 

The  electronic  degeneracy  in  a  molecule  forces  departures 
from  the  symmetry  on  which  the  degeneracy  was  evaluated  [51,52]. 
Therefore,  for  non-linear  molecules,  a  nuclear  arrangement  that 
will  lead  to  electronic  degeneracies  will  automatically  have  at 
least  one  non-totally  symmetric  vibrational  coordinate  that  will 
split  the  degenerate  state  and  destroy  the  symmetry.  When 
degeneracy  occurs,  only  those  systems  where  the  number  of 
electrons  either  completely  fill  or  half-fill  sets  of  degenerate 
one-electron  levels  will  not  Jahn-Teller  distort.  As  this  is  not 
the  case  of  either  Pd23  or  Pd4  structures,  the  nature  of  the 
distortion  should  be  evaluated. 

The  large  stability  of  the  Ijj  and  T^j  structures  for  the  13 
and  4  atom  clusters  respectively,  allowed  us  to  obtain  these 
structures  through  configuration  averaging.  This,  of  course,  is 
an  intermediate  stage  we  examine  to  make  comparisons  with  the 
work  of  others,  and  one  which  will  allow  us  to  examine  the  Jahn- 
Teller  distorsions  in  some  detail.  Of  interest  is  that  we  were 
unable  by  any  means  to  successfully  obtain  the  full  structure 
for  Rh23.  Rather,  a  near  ^5v  symmetry  is  achieved  from  both 
and  Ojj  starting  geometries  [25].  We  have  not  been  able  to  analyze 
quantitatively  the  steps  and  the  normal  modes  involved  in  the 
distortion  of  the  Rhj^3  molecule  from  to  Cg^,  although  we  could 
for  Rhi3'^  [25]. 

In  order  to  analyze  this  effect  for  Pd,  and  to  learn  about 
distortions  that  may  lower  the  energy  through  the  splitting  of 
the  degenerate  orbitals,  different  initial  geometries  have  been 


11 


compared,  which  are  related  through  a  normal  mode  to  the  higher 
symmetry  ones. 

a-  Icosahedral  Pd^3  clusters. 

According  to  the  electronic  distribution  depicted  in  Fig  3, 
the  Jahn-Teller  theorem  requires  a  splitting  of  the  (t2g)^  =  ^T2g 
+  %g  +  ^Ag  degeneracy,  leading  to  orbitals  of  e2g  and  a2g 
symmetry  in  Dg^.  Correspondingly,  the  ^T2g  state  reduces  to  ^E2g 
and  ^A2g,  the  ^Hg  state  to  states  of  ^A^g,  ^E^g  and  ^E2g  symmetry 
and  ^Ag  becomes  ^A^g,  see  Fig.  4.  The  relative  positions  of  the 
e2g  and  a2g  orbitals  in  energy  depend  on  the  final  geometry.  The 
simpler  distortions  that  can  be  idealized  imply  an  axial 
elongation  along  the  principal  axis  in  Dg^j  symmetry  (z  axis.  Fig 
2)  or  an  axial  compression  along  the  same  axis,  to  give  prolate 
or  oblate  ellipsoids  respectively,  from  an  initial  spheroid 
structure.  Depending  on  the  characteristics  of  the  distortion, 
the  e2g  orbital  will  be  above  or  below  the  a2g  for  the  elongated 
or  compressed  geometries  respectively. 

In  the  case  of  Pdj^3,  the  prolate  or  oblate  distorted 
starting  structures  returns  to  the  perfect  after  optimizing 
the  geometry.  Starting  structures  resulting  from  prolate  or 
oblate  distortions  along  a  3  fold  axis  in  03^1  symmetry,  or  a  2- 
fold  axis  in  02^  symmetry,  also  lead  to  the  perfect  Ijj. 

After  several  trials,  a  related  structure  was  found,  with  lower 
energy  than  the  initial  It  is  the  result  of  an  elongation 

of  the  Pd 2^3  Ijj  molecule  along  the  principal  axis  in  Dg^j  symmetry, 
coupled  to  a  compression  along  the  same  axis  that  shortens  the 
distance  between  oppossite  pentagons  (Fig.  5) .  This  corresponds 
to  an  hg  mode  in  Ijj  symmetry. 

The  molecular  orbital  distribution  in  the  Dg^j  structure 
shows  the  e2g  orbital  below  the  a2g  (Fig. 3).  The  lowest  energy 
corresponds  to  a  singlet  (^A^g) ,  followed  by  a  triplet  (^E2g)  , 
which  lies  0.16  eV  above.  The  next  states  are  a  triplet  and  a 
singlet  which  lie  0.31  and  0.49  eV,  respectively,  higher  in 
energy  than  the  GS,  see  Fig  5. 

The  distortion  of  the  geometry  can  be  described  as  a  first 
order  Jahn-Teller  effect  from  the  ^T2g,  in  Ij^  symmetry,  to  the 
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^E2g  in  05^,  through  an  hg  mode,  with  an  energy  decrease  of  0.9 
Kcal.  Further  decay  to  the  singlet  decreases  the  energy 

through  spin  pairing  to  give  a  closed-  shell  structure  (Fig.  5) . 
The  GS  of  the  Jahn-Teller  distorted  structure  is  4.68  Kcal  below 
the  spherical  ^T2g.  An  alternative  conceptual  description  may 
involve  the  lowest  singlets  in  and  05^  symmetries.  An  initial 
activation  from  the  triplet  to  the  singlet  implies,  however,  an 
energy  barrier  of  0.21  eV.  We  stress  conceptual  as  the  Jahn- 
Teller  effect  is  a  theoretical  construction  to  describe 
distortions,  not  a  physical  one.  One  does  not  observe 
distorting  to  Dg^,  one  simply  observes  Dg^. 

A  Mulliken  population  analysis  of  the  Dg^j  GS  structure  shows 
that  a  negative  charge  (-0.154  au)  is  still  located  on  the  atom 
in  the  center,  while  the  positive  charge  is  now  more  concentrated 
on  the  apical  atoms  (+0.149  au  on  each)  .  This  allows  one  to  infer 
a  larger  catalytic  activity  of  these  centers. 

Although  the  highly  symmetric  structures  are  characterized 
by  a  non-zero  magnetic  moment,  the  stabilization  through  a  Jahn- 
Teller  distortion,  followed  by  spin  pairing,  results  in  a 
prediction  of  diamagnetic  Pd23  clusters. 

The  consecutive  geometric  distortion-  followed  by  electron 
pairing  shows  the  importance  of  this  interplay;  the  geometry  of 
the  final  (distorted)  state,  which  is  a  consequence  of  the 
electronic  distribution  in  the  higher  symmetry  one,  leads  to  spin 
quenching . 

Providing  that  there  is  not  a  large  change  in  the  geometry 
and  that  the  molecular  orbital  distribution  is  nearly  kept,  the 
positively  charged  Pd23'‘’  structure  would  give  a  perfect 
icosahedron  which  would  not  distort,  due  to  the  equivalent 
filling  of  the  degenerate  t2g  HOMO  orbitals.  The  optimization  of 
the  geometry  gives  a  perfect  icosahedron  (Fig  2),  with  an 
interatomic  distance  of  2.63  A  and  a  molecular  orbital 
distribution  similar  to  that  in  the  neutral  but  for  the  t2g  -  hg 
inversion  and  the  larger  HOMO-LUMO  gap  (Fig  6)  .  The  GS  is  a 
quartet  ^Ag,  followed  by  a  doublet  ^Hg,  which  is  0.22  eV  higher 
in  energy. 
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Pdia"*^  clusters  would  be  predicted  to  have  non-  zero 
magnetic  moments. 

b-  Tetrahedral  Pd^  clusters. 

The  tetrahedral  Pd4  clusters  are  predicted  to  have  ^T2 
symmetry,  and  they  also  suffer  a  first  order  Jahn-Teller 
distortion. 

After  examining  several  structures,  we  find  that  the  most 
stable  structure,  depicted  in  Fig  7,  corresponds  to  a  ^A2  in  02^ 
symmetry,  and  is  the  result  of  a  Jahn-Teller  distortion  of  the 
^T2  in  T(j  symmetry  through  an  e  mode.  There  is  no  extra 
stabilization  through  spin  pairing.  The  GS  in  €2^  is  a  non¬ 
degenerate  triplet,  with  the  lowest  singlet,  ^A2,  0.38  eV  higher 
in  energy.  The  energy  associated  to  the  distortion  is  0.08  eV 
(1.84  Real) . 

The  Pd4  clusters  are,  thence,  paramagnetic. 

4-  CONCLUSIONS  - 

In  agreement  with  previous  results  on  Rhi3  and  Rh4  clusters, 
we  have  found  that  Icosahedral  and  Tetrahedral  geometries  are, 
respectively,  the  most  stable  ones  for  the  higher  symmetric 
structures.  Both  Ih  and  Td  structures  are  characterized  by  a 
large  concentration  of  triangular  faces,  which  are  the  minimun 
part  of  an  fee  (111)  structure,  associated  with  the  larger  packing 
density  and  the  lower  total  energy  for  this  crystalline  habit. 
This  is  accepted  as  a  valid  explanation  of  the  larger  stability 
of  these  geometries  for  the  cluster  sizes  analyzed.  The  final 
geometries  are,  however,  largely  determined  by  Jahn-Teller 
distortions,  that  we  predict  lower  the  and  T^j  geometries  to 
D5J  and  C2v  structures,  respectively. 

A  small  energy  difference  between  the  symmetrical  and  the 
distorted  structure  may  determine  the  spontaneous  evolution  to 
the  latter.  This  effect,  previously  found  in  Rh23  clusters,  does 
not  characterize  either  Pd23  or  Pd4  clusters,  and  the  different 
steps  involved  in  the  distortion  can  be  analyzed.  This  analysis, 
discussed  in  this  paper,  shows  the  close  relation  and  mutual 
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dependence  between  the  electronic  structure  and  the  geometry  of 
the  small  structures. 

The  electronic  configuration  of  the  final  structure,  which 
we  found  determined  by  a  Jahn-Teller  distortions  in  both  Pd4  and 
Pd3^3  defines  the  magnetic  moment  of  the  tiny  particles.  While  we 
calculate  that  paramagnetism  characterizes  Pd4  clusters,  the  most 
stable  Pd^j  structures  have  zero  magnetic  moment.  Recent  density 
functional  calculations  [50]  have  found  a  non-zero  magnetic 
moment  (-2jLtg)  ,  based  on  a  most  stable  Ijj  symmetry.  This  is  in 
agreement  with  our  results  for  this  structure.  But  our  final 
structure  is  a  distorted  cluster  of  geometry  and  has  a 
zero  magnetic  moment.  This  agrees  with  the  zero  magnetic  moment 
exper imenta 1 ly  f ound  [53]. 
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-  FIGURE  CAPTIONS  - 

-  Figure  1  -  The  frontier  molecular  orbitals  obtained  for  Pd4  in 
Tjj  symmetry  and  the  distorted  molecule  in  €2^  symmetry.  The 
distortion  that  relates  both  structures  is  indicated.  Energies 
are  in  Hartrees. 

-  Figure  2  -  Most  stable  Pd23  structure  as  result  from  a  full 
geometry  optimization  procedure,  starting  from  an  Ih  initial 
structure . The  symmetry  is  kept  when  4  electrons  are  averaged  in 
3  orbitals  in  the  CAHF  calculations.  More  details  are  given  in 
the  text. 

-  Figure  3  -  The  frontier  molecular  orbitals  obtained  for  Pd3^3  in 

symmetry  and  the  distorted  molecule  in  symmetry.  The 
distortion  that  relates  both  structures  is  indicated.  Energies 
are  in  Hartrees. 

-  Figure  4  -  Electronic  states  associated  to  the  ->  Dg^  Jahn- 
Teller  distortion  of  the  Pdj^g  structure. 

-  Figure  5  -  Distortion  of  the  Pd23  molecule,  in  Ij^  symmetry,  to 
the  most  stable  one  in  Dg^j  symmetry.  Geometrical  parameters  are 
indicated. 

-  Figure  6  -  The  frontier  molecular  orbitals  obtained  for  Pdj^g''' 
in  Ijj  symmetry,  from  CAHF  calculations  with  3  electrons  in  3 
orbitals.  Because  of  the  half  filling  of  degenerate  orbitals,  the 
molecule  does  not  Jahn-Teller  distort.  The  symmetry  of  the 
resulting  state  is  ^Ag. 

-  Figure  7  -  Distortion  of  the  Pd4  molecule,  in  T^j  symmetry,  to 
the  most  satble  one  in  0.2^  symmetry.  Geometric  parameters  are 
indicated. 
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Table  1-  Absolute  values  of  the  total  energies  [Hartrees] 
calculated  for  Pd4  tetrahedral  structures  after  MRCIS 
calculations.  Electronic  states  are  also  indicated.  Different  M 
are  compared  for  the  different  geometries  that  result  from  RHF 
(1st  row)  or  CAHF  calculations  with  different  number  of  open 
shells  included  in  the  average.  4/3,  4/6  and  6/9  stands  for  4  e" 
in  3  orbitals,  4  e“  in  6  orbitals,  and  6  e”  in  9  orbitals, 
respectively . 

In  the  last  row,  calculations  for  the  bulk  interatomic 
distance  are  reported.  The  Mulliken  populations  (M.P.)  given  are 
for  the  triplet,  see  text. 
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Table  2-  The  total  valence  energies  [hartrees]  calculated 
for  Pd;^3  structures  after  MRCIS  calculations.  Electronic  states 
are  also  indicated.  Different  M  are  compared  for  the  geometries 
associated  with  the  Ijj,  Oj^  and  D3j^  symmetries.  Details  on  the 
calculations  are  given  in  the  text.  The  Mulliken  populations 
(M.P.)  are  reported  for  the  central  atom  and  correspond  to  the 
most  stable  M  for  each  geometry. 


-0.21092 


02  -  -  0.24210 

^2 -  0.25206 

Bi  -  -  0.26038 


Ai  -H-  -  0.34588 

-0.3501  1  -44- 


-0.36348 


Bi  _4—  -  0.36568 
a1  -  0.36717 

B2  -  0.37226 


-0.41225  -44-  -44 - H-  ^2 


02 

Al 

01 


-  0.40590 

-  0.40819 

-  0.41038 


Pd4  ) 


2v 


3 


T2 


e 


3 


A2 


-0.26683 


eig 


Qlg 

^Z\} 

®1U 


hg 


a2g 


-0.29360 

-0.30094 


Qg 

tlu 


-0.41276  -JH- 

-0.42315  ^ 

-0.42546  -4V- 


+- 

+v 

-H- 


-0.44382 


^^13  i  U 


®2g 

eigV 


eiu 


-If- 


Pdl3  ,  D5J 


3 


0.21864 

0.21957 


0.24027 

0.24692 

0.25207 

0.26406 


0.41334 

0.42224 

0.42475 

0.42610 

0.42743 

0.42771 

0.43332 

0.44351 


AE[eV] 


.71 


S' 


-0,39590 


-0.42708 

-0.43422 


ag 

tiu 


-0.53781 

-0.54645 

-0.54698 


^ - 1-  t2g 

^  ^  _4_  h  , 


-  0.56617  -41-  -4i 


41 - H - 


Pd 


13 


Ih 


t 


